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Abstract 



A variant of the usual Lagrangian scheme is developed which describes both 
the equations of motion and the variational equations of a system. The re- 
quired (prolonged) Lagrangian is defined in an extended configuration space 
comprising both the original configurations of the system and all the virtual 
displacements joining any two integral curves. Our main result establishes 
that both the Euler-Lagrange equations and the corresponding variational 
equations of the original system can be viewed as the Lagrangian vector 
field associated with the first prolongation of the original LagrangianAfter 
discussing certain features of the formulation, we introduce the so-called in- 
herited constants of the motion and relate them to the Noether constants of 
the extended system. 
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I. INTRODUCTION 



The variational equations — this is Whittaker's [1] terminology — associated with dynam- 
ical systems are customarily obtained by linearizing the equations of motion around a par- 
ticular solution. Variational equations are important for understanding both stability and 
intcgrability issues [2-6] , and interesting since in general relativity and other metric theories 
of gravitation — as the JBD theory — they can be regarded as the equations of geodesic devi- 
ation, and also because they can be useful for describing linearized gravitation [7-10]. These 
equations are also useful for studying certain evolution equations [3,6,5,11,12]. Additionally, 
in chaos theory they are basic for defining the Liapunov spectra and the related Kolmogorov 
entropy [2,13]. 

This work deals with the variational equations associated to Lagrangian systems. This 
is not much of a limitation, however, as quantum mechanics, geometric control schemes, 
and field theories can all be described through a Lagrangian function. In all of these the 
problem of solving the variational equations is important [6,11,14]. We should also remember 
geodesies in Riemannian manifolds and their associated Jacobi fields whose properties have 
inspired much research and produced important results [15,16]. 

Though the equations of motion of the aforementioned theories stem from a Lagrangian 
through the Euler-Lagrange equations, the variational equations are normally formulated 
outside such scheme, the fact notwhitstanding that working within the Lagrangian formula- 
tion is convenient from a physical and a mathematical standpoint [2,5,17-19]. It is one of our 
aims here to describe a remedy to this situation, that is, we discuss a complete Lagrangian 
formulation of the variational equations and to take advantage of the description to get con- 
stants of motion. A feature of the approach is that the equations can be established making 
no reference whatsoever to any specific solution of the original equations, as is necessarily 
the case in the standard formulation. This formulation should also have importance in the 
theory of the Jacobi fields governing the transition from a geodesic to a nearby one in the 
calculus of variations — i.e. those vector fields which make the second variation to vanish 
identically excepting perhaps for boundary terms. The notion of Jacobi equations as an 
outcome of the second variation is in fact fairly more general than this and general formulae 
for the second variation and generalized Jacobi equations along critical sections have been 
already considered in the calculus of variations from a sort of structural, as has been called 
in [18], point of view (see also [20]). 

The paper is organized as follows. In section II we discuss the variational formulation 
first introduced in [2] . Section III deals with its the main features of the formulation, men- 
tioning some of its possible invariances. In section IV we exhibit how to every constant in the 
original problem it corresponds another one — what we call the inherited constant — valid in 
the variational system. In section V we prove Noether theorem for the variational equations, 
and establish under what conditions Noether constants can be regarded as inherited ones. 
We should pinpoint that the Noether theorem is able to reproduce the known constants as- 
sociated with symmetries of the original Lagrangian L, but that it also implies that to every 
n-parameter symmetry of L there additionally exist n new independent conserved quantities 
in the variational equations — which nevertheless can be trivial in some circumstances. In 
section VI, we use examples to pinpoint the kind of constants we may obtain by extending 
symmetries of the original Lagrangian [21]. One employs the simple case of ignorable coor- 
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dinates in classical mechanics to explain some further points on the application of Noether's 
theorem. The other uncovers a conserved vector that is valid in linearized gravitation in a 
vacuum. Section VI contains the conclusions and some final remarks. 



II. THE LAGRANGIAN FOR THE VARIATIONAL EQUATIONS 

Let us consider any system described in terms of a Lagrangian; that is, described using a 
real- valued function L defined in the tangent bundle TQ of the A^^-dimensional configuration 
space Q of the system, we additionally assume that L is non degenerate. The equations of 
motion of the system follow from 

The solutions of these Euler-Lagrange equations is the dynamics of the system. The formu- 
lation is described, excepting in section V and in one of the examples of section VI, under 
the implicit assumption that the parameter t is one-dimensional and it is usually referred to 
as the time; but it can actually be any other parameter, as the arc-lenght for geodesies in 
Riemannian manifolds, or even a finite set of parameters as in field theory. In such cases the 
necessary changes are easily done [16,23,21] as the reader can verify for herself (see section 
V and VI for specific examples). 

To describe deviations from the dynamics — the realm of the Jacobi equations — let us 
consider an augmented configuration space D comprising all the original configurations of the 
dynamical system plus all the elements of the Jacobi field that join two of its integral curves. 
The elements of D can hence be coordinatized by pairs (g", e") where the e" stand for the 
"virtual displacements" or deviations connecting two solutions of (1). In local coordinates 
we thus set q'"" = + and q'"- = + where is assumed "small". Formally we 
may regard D as the double tangent bundle of the configuration manifold of the original 
system composed with a canonical "fiip" mapping a reordering the local coordinates as 
{q,e,q,e) A {q,q,e,e) [22]. 

We can now define the new Lagrangian as a function defined on the tangent bundle TD 
of the D'Alambert configuration manifold and such that 

dL dL 

^) = + Q^ii^ = (2) 

where we have written the definition in TD^s local coordinates, we are using the summation 
convention and have taken the opportunity to define the operator V^. Thus 7 can be 
regarded as a directional derivative of L along a virtual configuration, or as the effect on 
L of an operator that "lifts it a little" out of the original space [24]. Formally 7 can be 
regarded as a prolongation [5,14] of L [2,17]. Note also that TD is the natural domain of 
definition of 7 since the virtual displacements are defined in the double tangent bundle of 
Q, and that the mapping a is well defined since the e's are elements of the tangent space 
TqQ at the point q in the configuration space Q. The A^^-dimensional object e hence plays 
the role of the variational field associated with the dynamical system's solutions [1,5]. We 
should pinpoint that despite the fact that we work in a specific chart of the configuration 
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manifold all the results can be translated into intrinsic language, as has been done in [5,22]. 
Be warned that the preliminary work reported in section 4 of [5] has been revised in [22]. 

Given the Lagrangian 7, we can immediately write down the 2N associated Euler- 
Lagrange equations of motion 



dt V ae« ) ae» 



and 



d f 97 
It \d^) ~ dq^ 



0. 



(3) 



(4) 



Using definition (2) in (3), we immediately obtain the Lagrange equations of the original 
system {i.e. we reproduce Eq. (1)) and, from (4), we get their associated linear variational 
equations [1,2,5,6] 



Mabi" + Cabi' + Kabe' = 0, a = 1, 2, . . . , TV, 



where the three objects M, C, K, are defined by the N x N matrices 



Mab^ 



a 



ab 



+ 



ab 



QqaQgb 

d_ (_^L_ 
dt ydq'^dq'' 
d_ f d^L 
dt \dq"-dq^J dq^'dq' 



_sn_ 

dq"-dq^ 



_dPL_ 
QqbQqa 

a,b = 



1,...,N. 



(5) 

(6) 
(7) 
(8) 



The variational equations (5) and the matrices M, C and K, are not explicitly time 
dependent unless L is it so from the start, a property not found in the standard approaches 
[1,6,11]. This happens because we are not linearizing around a particular solution but taking 
instead a less local point of view. Of course, at the end is all the same since for solving (3) 
we need to insert the solutions of (4). But the explicit time-independence of the coefficients 
in equations of motion (5) allows devising methods for analysing the variational equations 
analogous to those used in time-independent Lagrangian systems. Such lack of explicit 
f-dependence can be compared to what happens with the Lagrangian itself which is not 
regarded as time-dependent — unless it happens to be non autonomous — despite the fact 
that it depends on (g", g") and the solutions are always explicitly time-dependent. 



III. SOME FEATURES OF THE LAGRANGIAN FORMULATION. 

In this section we pinpoint the main features of the formulation; for a preliminary or a 
more mathematical outlook see, respectively, [2] and [22]. 

As it should be clear by now the function 7 is an entirely new Lagrangian which describes 
the original system plus its response to perturbations, i.e. what we have called the "virtual 
displacements" of the system. It is thus useful for studying both matters of stability and. 
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given the relationship of the variational equations with the Painleve test, also of integrability 

[1,2,5,3]. 

Moreover the Lagrangian formulation of the variational equations can be dressed in 
variational robes by using the function 7 to define the action functional 

^q{t),e{t)] = / j{q,q,e,e,t)dt (9) 

of the paths joining two (extended) configurations (q'i,ei) and (^'2,62) of the system at two 
given instants of time ti and t2. The statement of the variational principle is then just 
Hamilton's, that is, E should be extremal when the system follows its actual path in D 
[2,5]. In this way we can obtain Eq. (1) and Eq. (5) as a direct consequence of Hamilton's 
principle. The variational principle then becomes a compact formulation which may allow 
finding connections between the dynamical properties of the variational fields with other 
areas of physics or of mathematics [2,26,27]. 

As in any Lagrangian formulation, equations (1) and (5) are invariant under arbitrary 
point transformations — i.e. changes of coordinates — in Q, q"- = f"-{q,t), a — 1, . . . , N, 
where the /" are functions assumed to be invertible (det(rf/(g)) 7^ 0, g e Q). This result 
can be proved in the standard way [1,26]. 

It is well known [1,15,28] that the equations of motion (1) remain unchanged if we add 
to L the total time derivative of any function of q and t. This property is reflected in the 
invariance of Eqs. (3) and (4), obtained from 7, under the substitution 

f . ( ■ , f d'f{q,t) ., , d'f{q,t) \ , , df{q,t) .^ 

7(g, e, q, e) - ^{q, e, q, e) + [-g^q + ) e + -^e . (10) 

Note that the function / is arbitrary except for being required to depend only on q and on 
t. Of course, the equations of motion also remain unchanged when the total time derivative 
of any function, g, of e, q, and t is added to 7(5, e, q, e): 

, , . . .s dg(q,e.t) 
7(?,e,?,e) ^ 7(?,e,g,e) H (11) 

the proof of this property follows from the fact that the derivative in the right hand side 
of (11) satisfies identically Eqs. (3) and (4) [15,28]. Any transformation (10) is a particular 
case of (11) if we choose 

,i,,e,t)='J^e'. (12) 

This means that the extended system we are describing has a greater invariance than the 
original one; i.e. the extended system described by '-f{q,e,q,e,t) is invariant not just under 
the point transformations in Q explicitly mentioned above but also under the larger class of 
point transformations in D'Alambert's configuration space D [2,5,28]. 

It is to be noted that the e-derivatives of 7 reduce to corresponding ^-derivatives of L as 
follows 

|I = |£, and |I=|^1. (13) 

de"- dq°- de°- dq"- ^ ' 
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As a consequence of Eqs. (13), 7 is a first-order homogeneous function of the virtual dis- 
placements e" and velocities e" 

de"- ^ ^ 

and, therefore, it can be also written in the form 

where use have been made of the equations of motion (3), that is, relation (15) is valid along 
the integral curves associated with Eqs. (1) and (5). This equation also implies that the 
action S evaluated along such integral curves can be shown to be 

S=|^6» + Eo, (16) 

where Eq is the value of S at a certain convenient reference time to- 

There arc other relationships that can be written in terms of 7 and its derivatives. For 
example, the Hessian [16] of L, 

W = ^-^^e\' + -^e'^e" + ^-^^e^e'^ (17) 
can be expressed in terms of 7 as 

W = ^P,7 (18) 



5q"- dt \dq°- 



(19) 



where "D^ is the operator defined in (1), Sj/Sq is the functional derivative of 7 respect q 
[29,30]. We pinpoint that W can play an important role in variational problems since its 

sign distinguishes between minima, maxima and degenerate critical sections [16,18,28]. Note 
the similitude of the expression for W [Eq. (18)] with the definition of 7, this has interesting 
consequences for the theory of second variations when the action functionals are defined by 
first-order Lagrangians [18,19]. 



IV. CONSTANTS OF MOTION IN THE VARIATIONAL EQUATIONS. 

Let us consider any constant of motion, J{q, q), of the Lagrangian vector field Eq. (1), and 
let us evaluate it on two of the nearby solutions of Eq. (1), q(t) and q'(t) = q(t) + e{t). The 
difference between such constant quantities evaluated in nearby trajectories, j{q,q,e,e) = 
J{q', q') — J{q, q), is also trivially a constant, 

*(?^ = ^W,^')-ite<i)l = 0. (20) 
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This constant, j[q, e], which we call an inherited constant [5], can be also expressed as 



) 



(21) 



sort of directional derivative of J along a virtual configuration. Equation (21) tells 
us how, given both a solution of equations (1) and any one of its constants of motion, we can 
obtain a constant of motion for equations (5). This result establishes a direct relationship 
between integrals of motion in the variational equations, like j, to constants in the original 
equations (1). Related results are discussed in [6]. We need to pinpoint that, as in some 
non-linear evolution equations the integrals of motion J[q(t)] are functionals [29,30] — and 
not functions — of the solutions of (1), then the constant j has to be regarded also as a 
functional of the Jacobi fields e(t), given by [5,11] 



as has been shown in [6]; in Eq. (22) 5J[q{t)]/Sq{T) stands for the functional derivative of 
the functional J[q{t)]. 

At this point we recall the relationship between symmetries of L and the existence 
of integrals of motion in Lagrangian systems. This relation is summarized in Noether's 
theorem, which we, in the next section, proceed to extend to the existence of constants of 
motion in the variational equations. Such extension casts (21) in a different and perhaps 
more interesting form by relating it to a symmetry of the original Lagrangian L. But, as we 
mentioned before, the symmetries of the variational equations are larger than the symmetries 
of the original Lagrangian. Thus the Noether's theorem relates any continuous symmetry of 
7 — not necessarily coming directly from one of L — to a constant of motion in the variational 
equations. 



In the previous section we established the existence of what we called inherited constants 
of motion. As every isolating constant in a Lagrangian system come from a symmetry 
transformation [28,31], this has exhibited that to every symmetry of L there exist a related 
integral of motion in the variational equations. This also shows that any symmetry of L can 
also be regarded as a symmetry of 7; the converse of this statement is not necessarily true. 

In this section we formalize the just mentioned relation between constants of motion and 
symmetries of 7, making explicit the relation between the inherited constants of motion and 
the generators of the symmetry transformations of L. These results can be of particular 
importance in field-theoretic or in non-linear evolution equation applications [6,11,19,12]. 
Thinking on such field theoretic applications, in this section we assume that L depends on 
n parameters and not on a single one as we do in the rest of the paper. The n-parameters 
are further assumed to belong to a given n-dimensional differentiable manifold Vn as it would 
happen if they were points in the space-time continuum. 

Thus, let us suppose that the Lagrangian 7 is invariant under the r-parameter continuous 
group of transformations 




(22) 



V. SYMMETRIES AND INTEGRALS OF MOTION. 
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r^F^{t,q,e,w), r^FI{t,q,e,w), t^' ^ F-^{t,q,e,w), (23) 

in which the a;*, {s — 1, . . . , r) denote the r parameters of the group furthermore assumed 
independent of each other. The functions and Ft are parameterized by a;* G / (where 
/ is the set of parameters of the transformation group), and such that a;'' = 0, s = 1 . . .r, 
corresponds to the identity transformation, i.e. the group is continuously connected with 
the identity. The transformations (23) can be also written in the infinitesimal form 

r = ?" + C(^,?,e)c^^ (24) 
e" = e" + 77:(^,?,e)a;^ (25) 
i> = + ^)^s^ (26) 

where ^s{t,q,e), Csi^^ly^) r}g{t,q,e) are the infinitesimal generators of the symmetry 
transformation, given by 



That under the infinitesimal transformations (27-29) the system described by 7 remain 
invariant means that the action E [Eq. (9)] must behave under the symmetry transformations 
as [26] 

f „ dq" „ de"" -A , ^^t^'\ { n n ^e" ,A 
7 r , ^, 6^ ^, Ui;; = 7 ^, ^, . (30) 



The proof of the existence of the conserved quantities — details missing in the following 
proof can be filled in following the discussions in [26,31] — is direct. Begin with (30), differ- 
entiate it respect to uj^ noting that the right hand side is independent of such parameters, 
substitute a;* = after differentiation, and employ the equations of motion, to obtain 

^ = 0, (31) 
where the divergenceless tensor is given by 

fa = 7e + ^ (C - Q%Ca) + |r (^a - <^C) , (32) 

where we use the usual notation e^,y = de/dt". The result (32) is the main consequence 
of the Noether's theorem, establishing the connection between symmetries and conserved 
quantities at the level of the variational equations of a Lagrangian dynamical system. Let 
us point out that Eq. (32) was to be expected since 7 is indeed a Lagrangian. 

As the number of symmetries of the variational equations is larger than that of the 
original equations, we may ask under what conditions the Noether's constants (32) reduce 
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to the case of an inherited constant (21). If we regard the original constant of motion as 
coming also from a Noether symmetry, then 

J'a - Li^a + ^ (C - (33) 

the variational equations constant can always be written in terms of and of the operator 
"Dg as 

fa = T^eJ'a + ^ {t ' ^/.Cf ) " L ' ^ {C ' q'^sC!) ■ (34) 

where we have assumed the identity of the and of the ^'^ appearing in Eq. (33) with those 
in Eq. (34). It is quite clear then that not every constant can be regarded as an inherited 
constant. For this relation to be true two things are needed: i) that r]^ = i.e. that the 
transformation does not involve directly the e", and ii) that one or the other (or both) of 
the generators, Ca ^-iid Ca> be constants. 



VI. TWO EXAMPLES. 

In this section we give two examples of the use of Noether's theorem in the variational 
equations. 

A. Ignorable coordinates. 

Let us first consider a A^-degrees of freedom mechanical system described by the La- 
grangian, 

L{q, q) = rriabqaqb - V{q). (35) 

We assume this Lagrangian does not depend on the specific coordinate q^. The function 7 
is then 

/ 1 dniab . . dV\ . . 

^ ^ \ 2~dq~^''^^ ~ W)^^^ "^ab^a^fe- (36) 

and it neither depends on q"^. A one-parameter group of symmetries of both (35) and (36) 
is the translation 

g« = g«, ay^A, and g-^ = + (37) 

which is applied on Q, the original configuration space. The conserved quantity coming from 
such invariance is q^^s conjugate momentum pa = dL/dq^. We are now going to illustrate 
the relation between pa above [see Eq. (33)] and the constant in the variational equations 
[Eq. (34)]. The precise answer depends on how the transformation (37) is applied on D. 
As it is easy to convince oneself, the following two are the only independent possibilities 
allowed. 
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1. First, we apply (37) to and not to the "virtual displacements" or the time, i.e. 
= e"-, t = t, {(^ — 1, ("^ — 0, a ^ A; Tj"^ — 0; — 0). In this instance the conserved 
quantity is 

(1) _ ^7 _ druAb b b (oo^ 



and can be written as the inherited constant associated with pA, i.e. — V^pA- 

2. In the second case, transformation (37) is just applied to the e's and not to the q's 
or the t. We have = q'^; — + w, = e"', a ^ A; t — t {i.e. r]^ = 1, 
Tj"' = 0, a ^ A; = = 0). The conserved quantity is just the momentum conjugate 
to g^, j^A = d^/de^ = niAbq^ = Pa- 

It is worth emphasizing that this is the typical behaviour when we apply a symmetry of 
L to 7. Wc always obtain a new conserved quantity and a pre-existing (in L) one 
We pinpoint that it is not necessary that the constant j*^^) be inherited, as it happened in the 
present example; this might not even be possible. We must keep in mind that besides the 
inherited symmetries, 7 may have other symmetries of its own and therefore other associated 
constants. 



B. Conserved quantities in linearized gravitation in a vacuum. 

Let us investigate conserved quantities in a first-order theory of gravitation in a vacuum. 
The Lagrangian density [10] 

L = (39) 

with R°■^^^'^ the Riemann tensor, describes vacuum general relativity. The associated density 
is 

1 dR"'^^^ 

this Lagrangian describes linearized gravitation in a vacuum. 

It can be easily realized that under the following one parameter group of transformations 
of the spacetime 

^» = a;" + x^w, (41) 

the Lagrangian (39) changes as L = L/{1 + w)^. So, though it is not strictly invariant, L 
ends multiplied by a constant, thus there is a Noether tensor 9^ = —R"'^^'^Xa,fiXb,uX^/2 + 
R°'^^^Xa,^x^ which is not very interesting since it vanishes identically. Furthermore, the 
density 7 admits the following two (quasi) symmetry transformations: 
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For the first one, we just apply transformation (41), not changing the e's in any way. 
The Lagrangian 7 transforms as 7 = (1 + w)~^7, as a result the quantity 



= i^/^^e^..^-^ (42) 



is divergenceless [10]. 

2. The second invariance does not affect the spacetime coordinates x^ just transforms the 
"virtual" displacements e as 

e« = e» + e»w, (43) 

under this change the Lagrangian changes as 7= (l + u')~^7 and as a result the 
quantity j'* = d^/de'^^e"- = is also — trivially — divergenceless. This invariance does 
not give any useful new information. 

The first divergenceless tensor (42) , which do not stem from any constant in the original 
spacetime, is valid in linearized gravitation and most important remain constant irrespective 
of what spacetime we are linearizing around. 



VII. CONCLUSIONS 

The Lagrangian description discussed in this paper has led to the uncovering of various 
properties of the variational equations of Lagrangian systems, and allowed the application 
of Noether's theorem to them, thus relating the symmetries of a Lagrangian not only to its 
constants of motion but also to constants in its variational equations. These conclusions 
stem from the introduction of the Lagrangian 7 which, from a mathematical point of view, 
is the prolongation of L [5,22]. It is to be emphaziced that our results are vahd for field 
theories and non-linear evolution equations and can be also extended to the higher order 
Lagrangians which found application in relativistic theories [2,18]. Moreover, we can apply 
the techniques presented here in the framework of first-order variational principles on fibcrcd 
manifolds and their jet prolongations [19,35]. Our formulation may also possibly be used 
to study the so-called geoodular structure (a mathematical scheme that defines curvature 
in terms of properties of congruences of geodesies and their variational equations in affinely 
connected manifolds [36]). The extension of Noether's theorem to the variational equations 
presented here may also have some bearings in the study of the evolution of Stokes waves 
[12], and can be also used in geometric control theory [5,14]. 

On the other hand, the relevance of the variational equations for determining both sta- 
bility and integrability properties, makes this formulation important for the investigation 
of periodic orbits [3] and of solitonic solutions of nonlinear equations [12]. Our Lagrangian 
description of the Jacobi equations can be durthermore casted in Hamiltonian form. More- 
over, as the example B of section VI shows, the formulation lends itself to describe linearized 
gravitation. Let us mention that the mathematical setting of the concepts behind our for- 
mulation are explained more formally in [22] for the case of discrete Lagrangian mechanics. 
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